Momentum distribution of atoms in one-dimensional 
Raman cooling with arbitrary coherent population transfer 
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We derive an exact and analytical form for the cold-atom momentum distribution after a large 
number of one-dimensional (ID) Raman cooling cycles has been applied. Our result shows that one 
can select pulse profiles and lengths rather freely in order to obtain efficient cooling. Our approach 
takes into account optical pumping (resonant excitation followed by spontaneous emission) as an 
intergral part of the process, and it is noted that it does not affect the final momentum distribution. 
Thus subrecoil cooling is possible in general, and not just for carefully selected combinations of pulse 
forms and durations. 
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I. INTRODUCTION 

Cold atoms are a fundamental tool for studying degen- 
erate quantum gases [1], ultracold particle collisions [2, 
3], many-body effects [4], quantum computation [5, 6] 
and entanglement [7-9] , which then attracts attention to 
development of various cooling techniques. Laser cooling 
in combination with evaporative cooling [10] has made it 
possible to obtain Bose-Einstein condensation [11], and 
cooling of fermionic atoms via sympathetic cooling [12]. 
When quantum degeneracy is not required, laser cooling 
provides large numbers of cold atoms for a reasonable 
signal-to-noise ratio in atomic clocks [13, 14], as well as 
for unperturbed atomic transition frequencies in ultra- 
precise measurements [15]. 

Raman [16] and VSCPT [17] cooling have led to 
temperatures below the one-photon recoil limit in one- 
dimension (ID), but further experiments in 2D [18, 19] 
and 3D [18, 20] have certainly demonstrated the necessity 
for optimization. Also, proper design of state-insensitive 
traps is needed for quantum state engineering and preci- 
sion metrology [15, 21, 22]. In order to increase Raman- 
cooling efficiency, a tripod-type configuration of atomic 
levels was suggested for 2D cooling in Ref. [23, 24], as 
well as variants of coherent population transfer, based on 
Blackman [16, 18, 25], square [26, 27] or STIRAP [24, 28] 
pulses. A comparative analysis of subrecoil cooling is 
problematic, because it should predict the result of a 
large number of cooling cycles. On the other hand, even 
a qualitative analysis [26] given in the Levy-flight ap- 
proach [29] improved an earlier experiment [25] with a 
better fraction of atoms in the cold peak and a consider- 
ably simpler pulse sequence. In turn, our work reported 
in this paper provides a simple analytical solution for ID 
Raman cooling. 

Although single cycles in Raman cooling are usually 
well-defined, there is a large number of them, so that the 



cooling result becomes almost unpredictable. The coher- 
ent population transfer via Raman cycles, implemented 
already in many ways [16, 26, 28], can be obtained in 
numerous ways and it is not efficient to test them one 
by one. Thus the possibility to evaluate the cooling ef- 
ficiency quickly and in an simple manner is required. In 
addition, the efficiency of the cooling is affected by spon- 
taneous emission which involves a momentum broaden- 
ing of the cold atomic ensemble. This aspect must be 
accounted for when looking for efficient pulse properties. 
This paper is organized as follows. In Sec. II we de- 
fine the first step of a Raman cooling cycle with most 
general the transfer probability. The second step is de- 
scribed in Sec. Ill with an arbitrary angular distribution 
of spontaneous emission. We show that an analytical so- 
lution of the steady-state momentum distribution exists 
in general, first by deriving in Sec. IV an expression for 
it after a large number of cooling cycles have been ap- 
plied. The solution itself is derived in the case of most 
effective transfer profiles in Sec. V, and we conclude our 
presentation with a short discussion on the approach and 
its implications in Sec. VI. 



II. THE FIRST STEP OF COOLING CYCLE 

Figure 1 shows a three-level A-type atom that travels 
along the Oz axis, being prepared in the ground state 
|1). The momentum distribution of atoms is given by the 
probability amplitude a{p) of the atomic momentum p = 
Pz- The |1) state is coupled by a pump laser beam to the 
upper state |2), which in turn is coupled to the adjacent 
ground state |3) by a Stokes laser beam. The first step of 
a cooling cycle consists of a momentum-selective transfer 
accomplished through the two-photon transition |1) -H- 
|3). Defining the probability f{p) of the transfer, the 
number A7V of transferred atoms is given by 



A7V = /(p)|a(p)|2 



(1) 
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SO that the number of atoms left in the |1) state are 

\a,{p)\' = il^f{p))\a{p)\'. (2) 



Note that spontaneous decay from the upper level |2) is 
neglected here. Numerous studies of different pulse en- 
velopes such as Blackman [16, 18, 25], square [26, 27] and, 
recently, STIRAP [24, 28] pulses, have shown that the un- 
desirable contribution of spontaneous decay during typi- 
cal pulse durations is suppressed by the large frequency 
offsets of both laser beams from one-photon resonances. 
Thereby only the two-photon excitation is included and 
thus the width of resonant-momentum group is not lim- 
ited by the upper-state natural linewidth, making a deep 
subrecoil cooling feasible. For denser samples one needs 
to be concerned about collisional effects [2, 30, 31] and 
photoassociation resonances [3] , but these can be avoided 
e.g. by detuning the laser fields to the blue side of the 
transition [32, 33] as shown in Fig. 1. 

The number of lasers included into each laser beam 
as well as time-dependent changes in beam properties 
are not specified here. We are only interested in their 
contribution to the transfer momentum profile which is 
described by the function f{p). Independently whether 
a laser beam is provided by several lasers or not, let us 
assign it an overall wave vector: ki for the laser beam 
near-resonant with the |1) -f-)- |2) transition, and k2 for 
the one near-resonant with the |2) <-> |3) transition. In 
ganeral, Raman cooling uses one laser beam configura- 
tion to cool atoms with p < 0, and another one to cool 
atoms with p > 0. The first configuration is defined by 
ki = kie^, k2 = —^262 (shown in Fig. 1), and the sec- 
ond one by ki — — fciGz, k2 — k2ez- Depending on the 
configuration, atoms transferred to the |3) state gain a 
momentum shift of zLh{ki + k2), so that, in accordance 
with Eq. (1), the momentum distribution in state |3) is 
given by 

|a3(p)|' = f{p T Hk, + k2)MpTHki + fe))!^ (3) 

Here, signs "— " and "-I-" correspond to the first and sec- 
ond laser beam configuration, respectively. 



III. THE SECOND STEP OF COOLING CYCLE 




FIG. 1. A pump laser beam of wave vector ki couples tran- 
sition |1) <->■ 1 2) of a three-level A-type atom, while a Stokes 
laser beam of k2 couples the |2) -o |3) transition. A popula- 
tion transfer from the jl) to |3) ground state occurs through 
a two-photon transition. 



where N{u) is the angular distribution of spontaneously 
emitted photons. With the help of Eqs. (2) and (3), a 
change in the atomic population for the complete cooling 
cycle (steps 1 and 2 combined) is given by 

\a'ip)\'^{l-fip)Mp)\' 

+ / N{u)f(j)^hki+hkiu)\a{pT!T'ki +hkiu)\'^du. 

(5) 

Here, the left-hand side includes all the atoms of momen- 
tum p after the cooling cycle, whereas the right-hand side 
consists of two atomic parts: atoms remaining in the orig- 
inal state and those involved in the cooling process. A 
transferred atom starts from the |1) state at momentum 
Po = P -f Kki -\- hkiu, and subsequently experiences the 
total momentum shift 



6p = p-po = ±hki{lTu), 



(6) 



where the sign again depends on the laser beam configu- 
ration. 



The second step of a cooling cycle returns atoms into 
the |1) state by optical pumping. This is done by switch- 
ing the pump laser in Fig. 1 off and tuning the Stokes 
laser into resonance with the |2) ^ |3) transition. As 
an atom is excited from |3) to the upper state, its mo- 
mentum p' is shifted by hk2, becoming p' =p hk2, with 
the sign of change depending on the laser beam config- 
uration. Then the atom decays into the |1) state and 
a spontaneously emitted photon carries away a momen- 
tum kick of hkiu, modifying the atomic momentum as 
p = p' =F hk2 — hkiu. 

Consequently, the atomic population in state |1) after 
this step takes the form 



IV. LARGE NUMBER OF COOLING CYCLES 



|a'(p)|^ = |ai(p)|^ 



N{u)\a3{p± hk2 + hkiu)\^du, 

(4) 



To simplify equations, we set the scaling hki 
the nth cooling cycle is written as 



1. So 



Pn{p) - Pn-l{p) = -fn{p)Pn-l{p) 

+ N{u)fnipTl + u)pn-i{pTl + u)du, (7) 

where p„_i(p) and Pn{p) are the atomic populations 
|a(p)|^ before and after the nth cycle, respectively; fn{p) 
is the transfer probability. 

Equation (7) is of a special role in understanding the 
Raman cooling properties. As such a relationship is at- 
tained, the usual way to solve the issue is to subsequently 
specify transfer probabilities fn{p) a-nd angular distribu- 
tion N{u), and proceed with the investigation in the way 



of numerical calculations. Instead, we do not continue by 
specifying /„(p) and N{u), but consider repeated cooling 
cycles in general. Conditions when the atomic population 
Pn{p) tends to the steady-state solution p{p) as n — > cxd 
are discussed in Sec. V. 

To approach the steady-state solution p(p), wc consider 
a large number N of cooling cycles, following iVg initial 
cycles. The sum of Eq. (7) over N is given by 

No+N , 

PNo+n{p)~ PNaip)^ Yl i-fn{p)Pn-l{p) 
n=No + l ^ 

+ / N{u)fnipTl + u)pn-iipTl + u)duy (8) 

Once A'o is large enough, the atomic population ap- 
proaches the steady-state solution, i.e., Pnip) ~ p{p)j 
so that the left-hand side of Eq. (8) goes to zero, and 
equation (8) leads to inequality 

No+N , 1 

-/n(p)p(p)')«0. (9) 

In general, there are k different profiles for the first 
laser configuration, i^i(p), . . . ,Ffe(p), and the same 
number of profiles for the second laser configuration, 
Gi(p), . . . , Gfc(p). As a result, the sum (9) is split into 
2fc sums with similar transfer profiles, 

Y\\ N{u)Fi{p - 1 + u)p{p - 1 + u)du - Fi{p)p{p)\ , 

(10) 
and 



E 



N{u)Gi{p -Kl + u)p{p + 1 + u)du - Gi{p)p{p) 



(11) 

depending on the laser beam configuration. Since N is 
large enough, the number of cooling cycles in each of 2fc 
sums approximately equals N/{2k). After the substitu- 
tion of Eqs. (10) and (11) for Z = 1, . . . , /e into Eq. (9), 
the latter is written as 



N{u)F{p - 1 -t- u)p{p -1 + u)du - F{p)p{p) 



+ / N{u)G{p + 1 + u)p{p +l + u)du- G{p)p{p) = 0. 



(12) 

Here, Fiji) and G{p) are the average transfer probabilities 
for the first and the second laser configurations, respec- 
tively: 



F{p)^\J2Fi^P)^ G{p)^\J2Gi{p)- 



(13) 



The change from Eq. (9) to Eq. (12) is accompanied by 
change from an approximate equality to a strict one when 
A'^o and A^ go to infinity. 



V. THE MOST EFFICIENT EXCITATION 
PROFILES 

As follows from Eq. (6), a cooling cycle in the first laser 
beam configuration pushes the momentum of an atom 
along the Oz axis by a positive shift 5p. Thus transfer of 
atoms from the left-hand side of the momentum distri- 
bution causes cooling, and transfer from the right-hand 
side causes heating. To avoid the undesirable heating 
process, the transfer profile should preserve atoms with 
p > from the transfer, which is true for e.g. a Black- 
man and STIRAP-pulse envelope. Without focusing on 
a concrete momentum profile of the pulse, we impose a 
general condition 

F{p)>Q, ifp<0; F(p)=0, ifp>0, (14) 

in relation to the first laser beam configuration. For the 
second laser beam configuration the corresponding con- 
straint is written as 

G(p) = 0, ifp<0; G(p)>0, ifp>0. (15) 

In addition to efficient cooling, constraints (14) and 
(15) ensure that all atoms accumulate around the zero 
momentum. First of all, one can see that atoms within 
momentum region {—2hki,2hki) do not leave the region 
if conditions in Eqs. (14) and (15) are satisfied. Atoms 
from the rest of the ensemble are pushed closer and closer 
to the center of the distribution until they fall into this 
region (atoms on the left-hand side are pushed by means 
of the first laser configuration, and atoms on the right- 
hand side — by means of the second laser configuration) . 
Finally, the whole atomic ensemble will be transferred 
near the zero momentum, converging into a steady-state 
profile as the number of cooling cycles goes to infinity. 
Such a profile is described by the steady-state solution 
p{p)^ which in turn takes the form 



/9(p)=0, ifp<-2orp>2. 



(16) 



Whereas Eq. (16) only ensures that cooling exists, the 
cold-atom distribution is derived in Appendix A, where 
Eq. (A16) gives the corresponding steady-state solution: 



P{p) 



0, 

A 

A 

G{p) 
0, 



-i-p 



P<-2; 
N{u)du, pg(-2,0); 

I 

N{u)du, pe(G,2); 
p>2. 



(17) 



1=1 



1=1 



As follows from Eq. (17), the cold-atom distribution 
can be made arbitrary narrow once average transfer pro- 
files (13) have a proper form with a sharp dip at zero 
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FIG. 2. The cold-atom distribution against momentum p eval- 
uated for Raman cooling by STIRAP [28]. The analytical 
result (17) (solid line) shows the attainment of a cooling effi- 
ciency superior to that given in the reference (dashed line). 



momentum. A similar conclusion vifas pointed out in 
Ref. [28] for a population transfer produced by a STIRAP 
pulse. A comparison between the resulting cold-atom dis- 
tribution in Ref. [28] and that given by Eq. (17) is shown 
in Fig. 2. The analytical result (17) allows even more 
efficient cooling in comparison with Ref. [28] . 



VI. CONCLUSION 

We have presented the general result of ID Raman 
cooling as a simple analytical expression (17) which com- 
bines a large number TV of cooling cycles and gives the 
steady-state solution when N ^ oo. So far, Raman cool- 
ing has relied on the set of a specific order of profiles 
introduced by Kasevich and Chu in 1992 [16]; each pro- 
file in the set is narrower and closer to w = than the 
previous one in order to transfer most of the momen- 
tum distribution. Nevertheless, the resulting cold-atom 
distribution (17) is defined by averaged excitation pro- 
files (13) and is independent of the order of profiles. If 
needed, not only arbitrary order of excitation profiles can 
be applied with the same efficiency, but rather differ- 
ent profile envelops such as that produced by a STIRAP 
pulse [28], combining the excitation of a momentum- 
distribution wing with a narrow dip at w = 0. In addition, 
the steady-state solution is independent of any particular 
population transfer method, allowing one to observe the 
influence of spontaneous decay accompanying the second 
step in each cooling cycle. One can see that sponta- 
neous decay does not impose any constraint on neither 
the height or width of the cold peak, making deep sub- 
recoil cooling attainable. However, one should note that 
taking into account spontaneous decay through channel 
|2) o |3) gives more realistic situations, but leads to a 
significantly complicated consideration in the frame of 
the Wigner function [34], and we have not done it here. 
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Appendix A: Steady-state solution 

The steady-state solution p{p) in Eq. (12) is now mul- 
tiplied by either F{p) or G{p) transfer probability. Con- 
sequently, one can introduce functions 

yiip) ^ F{-p)p{-p), y2{p) = G{p)pip), (Al) 

which both are nonzero over range p € [0, 2) only, as 
follows from Eqs. (14), (15) and (16). Thus, functions 
2/1,2 (p) can be nonzero at p = due to the fact that p{p) 
can go to infinity as p — )■ 0. In terms of variables yi,2(p) 
Eq. (12) is written as 



N{u)yi{-p + 1 - u)du -~ 2/1 (-p) 



+ / N{u)y2{p + 1 + u)du - 2/2 (p) = 0. (A2) 

Although Eq. (A2) depends on two functions, 2/1 (p) 
and y2{p), it can be reduced for p E (0, 2) to an equation 
of one function 

v{p) = yi{2-p)+y2{p). (A3) 

First, we replace p by p — 2 in Eq. (A2): 



N{u)yi {-p + 3 - u)du - yi{2 - p) 



+ N{u)y2ip-l + u)du-y2{p-2) = 0, (A4) 

and then evaluate the sum of Eqs. (A2) and (A4): 
{2/1(2 -p)+ y2{p)} + {yi{-p) + 2/2(p - 2)} 

N{u){yi{-p+l - u) + y2ip + I + u)}du 



+ N{u){yi{-p + 3-u) + y2{p-l + u)}du. (A5) 

Excluding yi{~p) and 2/2 (p ~ 2), which both equal zero 
over p G (0, 2), the rest of Eq. (A5) is written in terms of 
function u(p) (A3) as 



v{p) 



i-p 



N{u)v{p+l + u)du 



+ N{u)v{p-l + u)du, (A6) 

n-p 



where we take into account that v{p) = if p ^ [0i2]. 

One can see that Eq. (A6) depends only on v(p) over 
range p G (0, 2). That is why if another function v{p) co- 
incides with v{p) for all p £ (0, 2), v{p) satisfies Eq. (A6) 
as well. At the same time, v{p) can take arbitrary values 
if p ^ (0,2). We assume that v{p) is a periodic function 
with period T = 2, hence Eq. (A6) applied for v(p) is 
written as 



x.)^/; 



N{u)v{p+ 1 + u)du. 



(A7) 



Here, we use relationship v(p + T) — w(p), and assume 
that v{nT) = yi(0) + 2/2(0) for all integer n. 

The periodic function v{p) is expanded into the Fourier 
series 



^(p) = E' 



^iknp 



Ofc. 



(A8) 



In terms of coefBcients a^ , Eq. ( A7) takes the form 
Xkak =0, Xk = l- (-1)^ f iV(u)e^'=""d7/. (A9) 

It is straightforward to verify that 

Xk=0, if fc = 0; Xk^O, iik^ 0; (AlO) 

hence, Eq. (A9) is equivalent to 

Uk = ASko, (All) 



where A is an arbitrary constant and Skn is Kronecker's 
delta. As a result, v{p) over range p G (0, 2) is given by 

v{p) ^ A. (A12) 

If p G (0,2), Eqs. (A2) and (A4) are written as 



2/2 (p) 



i-p 



N{u)v{p + 1 + u)du, (A13) 



yi{2-p)= N{u)v{p-l + u)du. (A14) 

Ji-p 

Hence, functions 2/1,2 (p) £^re given by 

yiip) ^ A / N{u)du, y2{p) = A / N{u)du. 
Jp-i J-i 

(A15) 

Finally, equations in Eq. (Al) give the steady-state solu- 
tion 



0, 



P<-2; 



p{p) = { 



A r^ 

-— / N{u)du, pe{-2,0); 
= } ^(^) J-^^-P (A16) 



A r^-p 



G{p) J., 
0, 



N{u)du, pG(0,2); 
P>2; 



where F{p), G{p) follow from Eq. (13). 
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